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Direct numerical simulations (DNS) of spatially developing turbulent boundary layers 
over riblets with a broad range of riblet spacings are conducted to investigate the effects of 
riblets on skin friction at high speeds. Zero-pressure gradient boundary layers under two 
flow conditions (Mach 2.5 with T w /T r = 1 and Mach 7.2 with T w /T r = 0.5) are considered. 
The DNS results show that the drag-reduction curve (A Cf/Cf vs if) at both supersonic 
speeds follows the trend of low-speed data and consists of a ‘viscous’ regime for small riblet 
size, a ‘breakdown’ regime with optimal drag reduction, and a ‘drag-increasing’ regime for 
larger riblet sizes. At if « 10 (corresponding to s + « 20 for the current triangular riblets), 
drag reduction of approximately 7% is achieved at both Mach numbers, and confirms the 
observations of the few existing experiments under supersonic conditions. The Mach- 
number dependence of the drag-reduction curve occurs for riblet sizes that are larger than 
the optimal size, with smaller slopes of A Cf/Cf for larger freestream Mach numbers. The 
Reynolds analogy holds with 2 Ch/Cf approximately equal to that of flat plates for both 
drag-reducing and drag-increasing configurations. 


Nomenclature 


A surface area, m 1 2 

Cf Skin friction coefficient, dimensionless 

Ch Stanton number, dimensionless 

C p heat capacity at constant pressure, J/(K-kg) 

C v heat capacity at constant volume, J/(K-kg) 

D Drag, N 

h riblet height, nr 

L Computational domain size, m 
M Mach number, dimensionless 
Pr Prandtl number, Pr = 0.71, dimensionless 
R radius of riblet tip 

Reg Reynolds number based on momentum thickness and freestream viscosity, Reg = Po °fff° e , dimensionless 
Re, 52 Reynolds number based on momentum thickness and wall viscosity, Re$ 2 = , dimensionless 

Re T Reynolds number based on shear velocity and wall viscosity, Re T = pw f" r5 , dimensionless 
s riblet spacing, nr 

T temperature, K 

T r recovery temperature, T r = T 00 ( 1 + 0.9 * K 

u streamwise velocity, m/s 

u T friction velocity, m/s 

v spanwise velocity, m/s 
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w wall- normal velocity, m/s 

z T viscous length, z T = u w /u T , nr 

a riblet ridge angle, degree 

7 specific heat ratio, 7 = C p /C v , dimensionless 

5 boundary layer thickness, nr 

<5* displacement thickness, m 

k thermal conductivity, k = pC p /Pr, W/(m-K) 

9 momentum thickness, m 

H dynamic viscosity, p = 1.458 x 10 T+110 4 , kg/(m-s) 

v kinematic viscosity, v = pj p, nr 2 -s 
p density, kg/m 3 

Subscripts 

f quantities related to flat plate or clean surface 

i inflow station for the domain of principal direct numerical simulation 
n surface normal direction 

r quantities related to riblets 

rms root mean sqaure 

w wall variables 

x streamwise direction of the right-hand Cartesian coordinate 

y spanwise direction of the right-hand Cartesian coordinate 

z wall-normal direction of the right-hand Cartesian coordinate 

00 freestream variables 

Superscripts 

+ inner wall units 

(•) averaged variables 

(•) / perturbation from averaged variable 


I. Introduction 

Riblets are micro-grooves aligned with the flow direction, and they are considered one of the premier 
techniques for reducing turbulent drag . 1 Laboratory experiments and numerical simulations have shown that 
riblets of different geometries can reduce skin-friction drag by up to 10 percent in turbulent boundary layers 
(see reviews by Walsh , 2 Choi , 3 and Garcia-Mayoral & Jimenez 4 ). Moreover, riblets have been successfully 
applied to full aerodynamic configurations . 5,6 For example, the flight tests of a commercial airplane (Airbus 
320) with riblets over 70 percent of its surface imply an overall 2 percent drag reduction, based on the fuel 
savings obtained . 7 

The performance of riblets in reducing drag depends on their size. Typically, the riblet dimensions are 
normalized by viscous length scale, z T , to account for the effect of Reynolds number. In the subsonic regime, 
both experiments and direct numerical simulations (DNS) show that a typical curve of drag reduction as a 
function of riblet spacing consists of a ‘viscous’ regime where the reduction in drag is proportional to the 
riblet size, a ‘breakdown’ regime where the linear dependence breaks down and a minimum drag is reached 
at some optimal spacing, and a ‘drag-increasing’ regime where riblets behave as the k-type roughness and 
induce a drag increase (Figure l ). 8 The viscous regime exists in the limit of very small riblets when the 
riblet is within the viscous sublayer, and the slope of the drag curve in the viscous regime, m s , can be well 
quantified by the Stokes theory . 9 Within the breakdown regime, riblet effect saturates, and a minimum 
drag is reached for spacings s+ pt rts 10-20. The optimal spacing s~ ( \ pt limits the maximum achievable drag 
reduction for a given riblet geometry. As the riblets get even larger, they become drag-increasing rather than 
drag-decreasing, and the slope of A t/t 0 versus s + curve within the drag- increasing regime determines the 
robustness of a particular riblet geometry in reducing drag. To design riblets with optimal peak performance 
and robustness, it is critical to identify flow parameters that could influence the extent of each regime and 
understand the mechanisms that control the transition among the different riblct-spacing regimes. Garcia- 
Mayoral and Jimenez 4 provided a recent review of available experiments and theoretical studies of the 
different drag-reduction regimes, as well as the physical mechanisms proposed in the literature both for 
the viscous regime and for its breakdown. They showed that the groove cross section, A g (or l g = ^yA g ), 
which accounts for both riblet spacing, s, and depth, h (Figure 2), provides a better collapse of the drag- 
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reduction curve for riblets of different geometries than the riblet spacing, s, with the optimal drag reduction 
being l^ opt ~ 10.7 ± 1. Garcia-Mayoral and Jimenez 3 4 ’ 10 also hypothesized that the viscous breakdown 
is caused by the appearance of long spanwise rollers below z + ss 20. The spanwise rollers are spanwise 
coherent structures that extend several inter-riblet spacings and that are developed from a two-dimensional 
Kelvin-Helmholtz-like instability of the mean streamwise flow. 

Most previous studies on the different drag-reduction regimes of riblets are limited to incompressible 
flows, and, typically, channel flows. Very limited riblet data at supersonic speeds are available. The few data 
include wind tunnel investigations at Mach 2.97 by Robinson, 11 Mach 1.5 by Gaudet, 12 Mach 1.6, 2.0, and 
2.5 by Coustols and Cousteix, 13 as well as flight tests at Mach 1.2-1. 6 by Zuniga et al. 14 These studies tested 
riblets with s + ss 9-15 and reported skin friction drag reduction of A Cf/Cf ss 4-15%. However, none of the 
above studies have provided a simultaneous measurement of skin-friction drag and near-wall turbulent flow 
field to allow for a complete characterization of the flow features in each regime and the mechanisms that 
control transition among regimes. Duan and Choudhari 15 conducted DNS of high-speed turbulent boundary 
layers to examine the effect of riblets on skin friction under two flow conditions (Mach 2.5 with T w /T r = 1 
and Mach 7.2 with T w /T r = 0.5). They considered symmetric V-groove riblets with two rib spacings 
(s+ ss 20 and s + ss 40). Unlike most of the previous DNS, 4, 16-18 which focused on incompressible channel 
flows, they examined riblets in external, spatially-developing flows without any assumption of streamwise 
periodicity. An identical inflow was used for simulations with different riblet spacings to provide a cleaner 
setting to compare drag. Their DNS results were consistent with the few existing experimental observations 
and showed that a drag reduction of approximately 7% can be achieved for riblets with proper spacing at 
high speeds. Their study also helped to establish a correspondence between riblets in an incompressible 
channel flow and those over a flat plate, especially under a supersonic and hypersonic free stream. However, 
with most of the work by Duan and Choudhari 15 focused on comparing the turbulence statistics and flow 
visualizations between drag-reducing and drag-increasing configurations, a complete characterization of flow 
features across different drag-reduction regimes for high-speed flows is still lacking. 

One of the primary objectives of the current study is to extend the work by Duan and Choudhari 15 to 
obtain a DNS dataset of supersonic external flows over riblet surfaces across all the riblet-spacing regimes 
(‘viscous’, ‘breakdown’, and ‘drag-increasing’), and provide detailed data for understanding drag-reduction 
mechanisms and developing reduced-order models for drag prediction. The effects of compressibility on 
the drag-reduction curve is explored and the similarities and differences in the drag reduction mechanisms 
between incompressible and high-speed boundary layers are investigated. 

The paper is structured as follows. The flow conditions and numerical methods are outlined in Section II. 
Section III presents results pertaining to the effects of riblets on skin friction at high speeds, including the 
drag-reduction curve and drag-reduction mechanisms, followed by a summary in Section IV. 

II. Flow Conditions and Numerical Methodology 

The boundary layer edge conditions and wall parameters for the DNS are given in Table 1, which provides 
freestream Mach number, density, and temperature (M^, p^, and T^, respectively) as well as the boundary 
layer thickness and various Reynolds numbers at the inlet of the domain. We conduct DNS for two spatially 
developing boundary layers in the high-speed regime. One is at Mach 2.5, which is relevant to supersonic 
transports, and the other is at Mach 7.2, which falls within the hypersonic regime. 

Table 1. Bounday layer edge and wall parameters for the DNS. 


Moo 

Poo (kg/m 3 ) 

Too(K) 

T w (K) 

T w /T r 

(5i(mm) 

Reg 

Re T 

Res2 

2.5 

0.1 

270.4 

568.0 

1.0 

4.00 

1,777.5 

340.4 

1,038.9 

7.2 

0.07 

65.8 

352.0 

0.5 

9.06 

6,536.3 

392.4 

1,372.3 


The compressible Navier-Stokes equations are solved in generalized curvilinear coordinates. The working 
fluid is assumed to be an ideal gas with a linear (i.e. , Newtonian) stress-strain relation. The Fourier law is 
used to compute the heat flux terms. A 7th order weighted essentially non-oscillatory (WENO) scheme 19 
is used to compute the convective flux terms. This particular WENO algorithm combines a high order 
of accuracy with low dissipation and high bandwidth, which is desirable for time accurate simulations of 
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compressible turbulence. The resolution properties of this scheme are documented in several references. 20,21 
For the viscous flux terms, a 4th order central difference scheme is used. The 3rd order low storage Runge- 
Kutta scheme by Williamson 22 is used for time integration. 

To enable the simulation of a spatially evolving boundary layer, we conduct an independent auxiliary DNS 
of a spatially developing boundary layer which provides the inflow boundary condition for the principal DNS 
over the domain with riblets (Figure 3). The single auxiliary DNS is used to provide an identical inflow for the 
downstream simulations with and without riblets, thus isolating the influence of riblets on skin friction drag. 
To minimize the extent of upstream/downstream influence due to the surface variation, the riblets begin 
at a distance downstream of the inlet of the principal DNS and the grids consist of a smooth streamwise 
transition from the clean surface to the riblet surface (Figure 4a). 15 On the wall, no-slip conditions are 
applied for the three velocity components and an isothermal condition is used for the temperature, with the 
wall temperature listed in Table f. The wall temperature is prescribed to be nearly the recovery temperature, 
T r . for the Mach 2.5 cases and 0.5T r for the Mach 7.2 cases. At the top and outlet boundaries, unsteady 
non-reflecting boundary conditions based on Thompson 23 are imposed. Periodic boundary conditions are 
used in the spanwise direction. 

For each freestream condition, simulations with different surface topologies have been undertaken with 
riblet spacing the primary parameter as outlined in Table 2. To our knowledge, this is the first time that 
such a systematic parameter sweep has been undertaken for numerical riblets under supersonic speeds. For 
each riblet simulation, triangular riblets with ridge angle of 45° are used. To improve smoothness of the 
grids and avoid singularities when evaluating the Jacobian, the riblet tips and valleys are slightly rounded 
with with radius R/s ss 4%. The extent of rounding is within a range where no significant degradation in 
riblet performances is expected. 24,25 More importantly, as discussed by Garcia-Mayoral and Jimenez, 10 such 
a tip-rounding should not significantly affect the viscous breakdown, because they barely modify the groove 
cross section. A close-up view of riblet geometry and near-wall grids has been sketched in Figure 4b. 

For DNS cases with riblets, the riblets begin at a distance x ss 4.8 Si downstream of the inlet of the 
principal DNS, while upstream of this position the surface remains clean to guarantee a nearly unperturbed 
flow at the inlet. Moreover, the grids consist of a smooth streamwise transition from the clean surface to 
the riblet surface to minimize the extent of upstream/downstream influence due to the surface variation. 
The streamwise length for the transitional zone is approximately four times that of the riblet height, h, 
with a finer spacing across the transition length, as is shown in Figure 4a. The streamwise grid is uniform 
on both sides of the transitional zone. The riblet zone has a streamwise extent of more than 10 Si. The 
spanwise grids are clustered near the riblet tips and coarsened near riblet valleys, with 32 grid points on the 
surface of each riblet. The spanwise extent of the domain is approximately 400 wall units. Such a spanwise 
extent is significantly larger than the values used by Choi and Moin 16 and comparable with those used by 
Garcia-Mayoral and Jimenez. 4 


III. Results 

In this section, we present DNS results based on the cases listed in Table 2. The effects of riblets on 
skin-friction drag and other turbulence statistics are examined, and the comparison of flow physics among 
different Mach-number cases is conducted. Statistics are collected near a streamwise location of x a « 18 S t 
for the Mach 2.5 cases and x a ~ 15 for the Mach 7.2 cases, which is significantly downstream of the beginning 
of the riblet surface to avoid any significant effects due to surface transition. To obtain statistics, averages 
are calculated over time, the same spanwise locations over different riblets, and a streamwise window of 
[x a — 0.8 5i,x a + 0.8 <5j]. The time average spans more than 150 Si/Uoo- We use an overbar to indicate an 
averaged quantity and a prime to indicate perturbations from this average. Averages are calculated over 
time, the same spanwise locations over different riblets, and the streamwise window. Reducing the number 
of flow fields for the ensemble averaging by one third or increasing the length of streamwise window by 
eight times result in less than 0.5% difference in viscous drag. The origin of the 2 coordinate corresponds to 
the wall location for clean surfaces. For cases with riblets, the riblets are mounted on top of the otherwise 
clean surface with 2 = 0 corresponding to the riblet valley. This computational setup mimics the practical 
applications wherein riblets are molded into a lightweight plastic film with an adhesive backing and pressed 
into place on the surface of a vehicle. The influence of the surface transition from clean surface to riblet 
surface has been assessed by varying the streamwise transitional length, and less than 0.5% difference in 
drag is observed at the downstream location where statistics are collected. 
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Table 2. Simulation parameters for the direct numerical simulations of turbulent flows over 
riblets. The value of Ax corresponds to that in the uniform region (outside of the transitional 
zone). 


Case 

s+ 

l t 

N riblet 

L x /5i 

Ly/Si 

L z /Si 

Ax + 

+ 

<1 

Az min 

N x 

X Ny X 

N z 

M25Clean 

NA 

NA 

NA 

20.4 

1.6 

11.3 

10.3 

1.34 

0.28 

528 x 320 x 120 

M25sl0 

8.9 

4.5 

40 

20.4 

1.6 

11.3 

10.3 

0.13-0.51 

0.20-0.27 

528 x 1280 x 120 

M25s20 

18.2 

9.1 

20 

20.4 

1.6 

11.3 

10.3 

0.27-1.02 

0.20-0.27 

528 

X 

640 x 

120 

M25s30 

27.6 

13.8 

14 

20.4 

1.6 

11.3 

10.3 

0.39-1.53 

0.20-0.27 

528 

X 

448 x 

120 

M25s40 

37.3 

18.6 

10 

20.4 

1.6 

11.3 

10.3 

0.52-2.04 

0.13-0.28 

528 

X 

320 x 

120 

M25s50 

47.6 

23.8 

8 

20.4 

1.6 

11.3 

10.3 

0.65-2.55 

0.13-0.28 

528 

X 

256 x 

120 

M25s60 

58.2 

29.1 

7 

20.4 

1.6 

11.3 

10.3 

0.78-3.06 

0.13-0.28 

528 

X 

224 x 

120 

M72Clean 

NA 

NA 

NA 

16.8 

1.1 

14.5 

5.8 

1.18 

0.28 

960 

X 

320 x 

120 

M72s20 

20.8 

10.4 

20 

16.8 

1.1 

14.5 

5.8 

0.16-1.29 

0.19-0.52 

960 

X 

640 x 

120 

M72s30 

32.1 

16.0 

14 

16.8 

1.1 

14.5 

5.8 

0.24-1.93 

0.19-0.57 

960 

X 

448 x 

120 

M72s40 

43.8 

21.9 

10 

16.8 

1.1 

14.5 

5.8 

0.32-2.64 

0.19-0.62 

960 

X 

320 x 

120 

M72s50 

56.5 

28.2 

8 

16.8 

1.1 

14.5 

5.8 

0.45-2.95 

0.18-0.62 

960 

X 

256 x 

120 

M72s60 

59.9 

34.9 

7 

16.8 

1.1 

14.5 

5.8 

0.64-3.86 

0.18-0.62 

960 

X 

224 x 

120 


A. Drag Measurement 


The total skin-friction drag, D 1 is obtained as the integral of the product of the wall shear stress and the 
differential wetted area 

D= L/ w fe) ,,, dAr = <,Af = fiw fe) f Af (i) 

where A r denotes the wetted area of a riblet surface and Af = A r cos(a) = L x x L y is its corresponding 
clean-surface area (see Figure 5). 

The skin friction coefficient, Cf , is defined as 


C f = 


2PooUoo 


( 2 ) 


Figure 6 shows typical time history of the plane-averaged skin- friction coefficient at x/6i = 18 after the 
initial transient for the Mach 2.5 cases. Consistent with the results at low speeds (Figure 1), the performance 
of riblets is a function of riblet size and there exist a drag-reducing regime and a drag-increasing regime. 

Figure 7 further plots the drag-reduction curve, ACf/Cf (time and space-averaged skin-friction variation 
relative to the clean surface), as a function of riblet size, /+ = (A+) 1 / 2 , for the various freestream Mach 
numbers. For reference, the envelopes of low-speed experimental results compiled by Garcia-Mayoral and 
Jimenez 4 are included. Note that Z+ rather than s + is used to characterize the riblet size following the 
definition introduced by Garcia-Mayoral and Jimenez. 4, 10 Such a choice will facilitate the comparison with 
the low-speed data, since /+ has been shown to result in a better collapse of data at low speeds for riblets 
of different geometries. 4,10 Similar to the low-speed counterpart, the drag reduction for both supersonic 
Mach numbers attains its maximum at 1+ ss 10, with a maximum drag reduction of approximately 7%. 
The achieved drag reduction at both supersonic Mach numbers confirms the observations of the few existing 
experiments under supersonic conditions. 12-14 The similarity in the amount of drag reduction at ss 10 
among various Mach number cases may indicate that the dynamics of riblet flows in the breakdown regime 
is not significantly modified by compressibility effects up to a Mach number of 7.2. 

The Mach-number dependence of the drag-reduction curve seems to occur for riblet sizes that are larger 
than the optimal size l^ op t ~ 10. As the freestream Mach number increases, the slope of A Cf/Cf at 
> l+ opt , which determines the robustness of a particular riblet geometry in reducing drag, reduces. As a 
result, the optimal drag reduction deteriorates at a lower rate and the drag-reducing regime persists up to 
a larger riblet size for the higher Mach number case. 
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It should be noted that the above observations are only made based on the current data set. Additional 
data points in the vicinity of the maximum-drag-reduction location are necessary to confirm that l^ opt stays 
approximately around 10, especially for the Mach 7.2 case, for which there is no data between 0 < Z+ < 10.4 
within the current data set. Also, the extent of the ‘viscous’ regime and the exact location where the drag- 
reduction curve begins to deviate from a linear relation cannot be fully characterized based on the current 
DNS datasets alone and needs to be determined via additional DNS together with Stokes calculations. 9 

B. Near-wall Turbulence Structures and Drag Reduction Mechanisms 

In this subsection, near-wall turbulence structures are investigated in order to identify the similarities and 
differences in the mechanism by which riblets reduce drag between incompressible and supersonic flows. 

In a recent DNS study of riblets for incompressible flows, Garcia-Mayoral and Jimenez 4 found that 
the viscous breakdown is caused by the appearance of long spanwise vortical structures (referred to as 
spanwise rollers following Garcia-Mayoral and Jimenez) below (z — z t i p ) + ss 20 with typical streamwise 
wavelengths A+ ss 150 that are developed from a two-dimensional Kelvin-Helmholtz-like instability of the 
mean streamwise flow. Figure 8a and Figure 8b plot the instantaneous wall-normal velocity at (z — Zu p ) ~ 3 
to visualize the near-wall spanwise rollers for both supersonic Mach-number cases at l + ~ 10. The spanwise 
rollers that are perpendicular to the riblets appear for both Mach-number cases (e.g., at x + ss 1,200 for 
the Mach 2.5 case), although it is much more obvious for the Mach 7.2 case. The existence of the spanwise 
coherent structures may indicate a connection between the spanwise rollers and the viscous breakdown at 
high speeds, as in the low-speed flows. 

In terms of the drag-increase regime (corresponding to Z+ > Z+ opt ), Choi and Moin 16 argued that the 
deterioration of drag reduction for large riblet sizes is caused by the lodging (penetration) of the streamwise 
vortices inside the grooves so that a larger wetted area is exposed to the high-speed downwash induced 
by the streamwise vortices, resulting in a larger region of high shear rate. There exist similar streamwise 
vortices for high-speed flows, as visualized in Figure 9, and the characteristics of these vortices vary with 
the freestream Mach number and wall temperature. 26,2 ' To estimate the average size and location of the 
near- wall streamwise vortices, Figure 10 plots the root-mean-square (RMS) of streamwise vorticity, u) xrms , 
as a function of the wall-normal distance for Cases M25Clean and M72Clean. For both supersonic Mach 
numbers, u xrms shows a local minimum at z^ in ~ 5 and attains a local maximum at z+ oa , ~ 20, with 
z max = 1^.4 and z+ lin = 4.5 for Mach 2.5 and z^ ax = 23.2 and z+ nn = 4.9 for Mach 7.2. According to 
the streamwise vortex model proposed by Kim et al. 28 (Figure 11), the local maximum corresponds to the 
average location of the center of the streamwise vortex, and the local minimum corresponds to the average 
location of the edge of the vortex. The larger values of z+ aa , and (z+ aa , — z+ in ) for the Mach 7.2 case 
indicate that the streamwise vortices for the Mach 7.2 cases are of larger size and located farther from the 
wall compared with the Mach 2.5 case. As a result, for a given groove size, it is less likely for these vortices 
to lodge within or penetrate into the grooves, as shown in Figure 12. This may explain why the drag increase 
is less rapid for the Mach 7.2 case than the Mach 2.5 case for Z+ > l popt - 

C. Surface Heat Flux 

The design of effective vehicle surfaces that minimizes surface heat flux is of paramount importance for 
preventing overheating of hypersonic vehicles. Thus, one of the most important factors that determine the 
applicability of riblets to hypersonic vehicles is their influence on surface heat flux. Although heat-transfer 
management using riblets has been studied extensively for low-speed flows, there is some controversy about 
their effectiveness, 29 31 and very few studies exist for high-speed flows. In this subsection, the effect of riblets 
on surface heat flux is investigated for the Mach 7.2 case, which has a wall temperature of T w /T r = 0.5. 

Similar to the total skin- friction drag, the total surface heat flux, Q, over a riblet surface, A r , can be 
obtained as 



The Stanton number Ch are defined as 


C h 


— * 
Qw 


Pc<)^ooCp(T r T w ) 


( 4 ) 
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and the Reynolds analogy factor RA is defined as 


RA = 


2 C h 
C f 


( 5 ) 


Figure 13a shows the variation in Stanton number, Ch, as a function of riblet size, Z+, for the Mach 7.2 
case. Similar to the skin-friction drag, the surface heat flux reduces for small riblet sizes and increases for 
large riblets, with the maximum heat reduction at Z+ ss 10. The fact that a drag-reducing riblet surface also 
reduces surface heat flux is a feature that is desired in hypersonic applications, even if the reduction in heat 
transfer is modest. 

For drag-reducing and drag-increasing riblet surfaces, the variation in heat transfer is larger than that 
in viscous drag with the Reynolds analogy factor 1.16 < RA < 1.19 (Figure 13b). The Reynolds analogy 
factor for riblets assumes similar values as those for hypersonic flat-plate boundary layers. For comparison, 
the experimental data reviewed by Hopkins and Inouye 32 give 1.1 < RA <1.3 for boundary layers with 
T w /T r > 0.3. RA ss 1 indicates the validity of Reynolds analogy for riblet surfaces, at least for those with 
Z+ < 30. 

The relatively larger variation in heat transfer than that in viscous drag and the validity of Reynolds 
analogy across a wide range of riblet sizes is consistent with the findings by Benhalilou and Kasagi 33 and 
Stalio and Nobile 31 for low-speed boundary layers. However, given that the temperature is only a passive 
scalar for low-speed flows but is tightly coupled with the velocity field for high-speed flows, it may be difficult 
to draw any direct analogy between the two speed regimes. 


IV. Summary 

This paper presented direct numerical simulations of turbulent boundary layers over triangular riblets 
with a ridge angle of 45° at Mach 2.5 and Mach 7.2. The DNS results show that the drag-reduction 
curve (ACf/Cf vs riblet size) at both supersonic speeds follows the trend of low-speed data and consists 
of a ‘viscous’ regime for small riblet size, a ‘breakdown’ regime with optimal drag reduction, and a ‘drag- 
increasing’ regime for larger riblet sizes. At Z 4 " ss 10 (corresponding to s + ss 20 for the current triangular 
riblets), drag reduction of approximately 7% is achieved at both Mach numbers, and confirm the observations 
of the few existing experiments under supersonic conditions. The drag-reduction curve as a function of riblet 
spacing displays a noticeable dependence on the flow Mach number, with a weaker rate of reduction in drag 
benefit at the higher Mach number. 

In terms of drag-reduction mechanisms, the DNS data show the existence of spanwise coherent vortical 
structures at Z + « 10, similar to the finding by Garcia-Mayoral and Jimenez. 4 The reduction in the slope 
of A Cf/Cf for Z+ > l^_ opt at the larger Mach number appears to be related to the location of streamwise 
vortical structures relative to the riblet grooves. 

At Mach 7.2 with T w /T r = 0.5, it is found that the Reynolds analogy holds with 1.16 < RA = 2 Ch/Cf < 
1.19 for both drag-reducing and drag-increasing riblet surfaces. 

Future work will target additional computations to determine the optimum riblet spacing at the higher 
Mach number as well as to provide a better characterization of the extent of the ‘viscous’ (i.e. , linear) portion 
of the drag reduction curve at smaller riblet spacings. 
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Figure 1. Definition of the drag-reduction regimes observed over triangular riblets with 60° tip 
angle, as a function of the peak-to-peak distance s + . Adapted from Bechert et al. 8 
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Figure 2. Sketch of groove cross section A g . Adapted from Garcia-Mayoral and Jimenez. 4 
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Figure 3. Sketch of DNS domain for spatial simulations. 
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(a) Side view of the transition from clean surface to (b) Close-up view of riblet configuration 

riblets 


Figure 4. Computational mesh near riblets for M25s20. 



Figure 5. Sketch of riblet geometry. Adapted from Robert. 34 



tu / 5 , 


Figure 6. Time history of the skin-friction coefficient for various Mach 2.5 cases. 
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Figure 7. Drag-reduction results from DNS of supersonic flat-plate boundary layers with tri- 
angular riblets. The symbols represent DNS results with the error bars estimated based on the 
variability in drag reduction as a function of the size of the averaging window in x and t. The solid 
lines represent the envelope of the experimental data at low speeds compiled by Garcia- Mayoral 
and Jimenez. 4 



X + X + 

(a) M25s20 (b) M72s20 


Figure 8. Instantaneous wall-normal velocity at ( z — zu P ) + ~ 3 within the riblet zone of the 
simulation domain. The origin in x has been shifted arbitrarily. 



x + x + 

(a) M25Flat (b) M72Flat 


Figure 9. Instantaneous streamwise velocity at z + ~ 15 for M25Flat and M72Flat. The origin in 
x has been shifted arbitrarily. 
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Figure 10. The RMS of streamwise vorticity as a function of the distance away from wall for 
M25Clean and M72Clean. 



Figure 11. A near-wall streamwise vortex model (Adapted and modified from Kim et al. 28 ). In 
the figure, Zmax and z m in correspond to the local minimum and local maximum of the rms of 
streamwise vorticity shown in Figure 10. 



Figure 12. The RMS of streamwise vorticity as a function of wall-normal distance above the riblet 
valley for several riblet cases. 
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Figure 13. Stanton number and Reynolds analogy factor as a function of riblet size l + for the 
Mach 7.2 case with T w /T r = 0.5. 
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